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Abstract 

We present a nonperturbative approach which allows to evaluate the weak localization correction to the conductivity 
of disordered conductors in the presence of interactions. The effect of the electron-electron interaction on the 
magnetoconductance is described by the function A{t) exp{—f{t)). The dephasing time is determined only by f{t), 
and this time remains finite down to T = due to the electron-electron interactions. In order to establish the 
relation between our nonperturbative analysis and the perturbative results the effect of interaction on the pre- 
exponent A{t) should be taken into account. The dephasing time cannot be unambiguously determined from a 
perturbative calculation. 
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1. Introduction 

Recent experiments by Mohanty, Jariwala and 
Webb [1] attracted a lot of attention to the fun- 
damental role of interactions in disordered meso- 
scopic systems. These experiments strongly indi- 
cate that the low temperature saturation of the 
electron decoherence time [2,3] in disordered 
conductors has an intrinsic nature. Further analy- 
sis [4] allowed to rule out various experimental ar- 
tifacts and several theoretical explanations based 
on extrinsic effects. 

Can one expect the effect [1] theoretically? In [5] 
we developed a nonperturbative theoretical anal- 
ysis of the above phenomenon. We demonstrated 
that the low temperature saturation of the deco- 



herence rate in disordered mesoscopic sys- 

tems can have an intrinsic nature: it can be caused 
by the electron-electron interaction. The results [5] 
are in a good agreement with experimental find- 
ings [1]. We also argued [6] that this interaction- 
induced decoherence at low T has the same physi- 
cal nature as in the case of a quantum particle in- 
teracting with a bath of harmonic oscillators [7-9] . 

In contrast, in [10,11] and in some other papers it 
was argued that interaction-induced dephasing of 
electrons at T — > is not possible. Let us review the 
main arguments in favor of this conclusion [10,11]. 

One of such arguments is not necessarily related 
to electrons in a disordered metal. One can ar- 
gue that a particle with energy ~ T cannot excite 
harmonic oscillators with frequencies u > T and, 
hence, the latter can only provide some renormal- 
ization effects. Within this scenario electron scat- 
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tering on such a static potential is purely elastic 

and cannot lead to dcphasing. 

It is easy to observe that the above argument 
explicitly contradicts to the exact results obtained 
e.g. within the Caldeira-Leggett model [7]. In this 
model even at T = the off-diagonal elements of 
the particle density matrix decay at a finite length 
set by interaction with an effective environment. 
This effect is due to all high frequency modes, i.e. 
the picture is by no means "static". Our results 
[5] demonstrate that also in disordered conductors 
the high frequency "quantum" modes with u) > T 
do contribute to electron dcphasing. 

One can also argue that the electronic system 
can behave differently from a bosonic one [7] be- 
cause of the Pauli principle which restricts the elec- 
tron ability to exchange energy at low T. Again, 
this argument contradicts to the well known re- 
sults obtained for fermionic systems. E.g. tunnel- 
ing electrons exchange energy with the effective en- 
vironment (formed by other electrons in the leads) 
even at T = [13,14]. This so-called "P(E)-theory" 
[13] yields measurable consequences and was veri- 
fied in many experiments [15]. A close formal and 
physical similarity between the P(E)-theory and 
our analysis [5] is discussed in [16]. 

A more formal argument is based on a pertur- 
bative calculation [11]. The authors [11] argued 
that their results explicitly contradict to ours and, 
hence, the latter cannot be correct. A direct com- 
parison between the above results is not quite sim- 
ple. The main reason for that is that our calculation 
[5] is essentially nonperturbative. It was performed 
with the exponential accuracy which is sufficient 
to determine and the weak localization correc- 
tion to the conductance da in the leading approxi- 
mation for all values of the magnetic field. Aleiner 
et al. [11] considered the limit of strong magnetic 
fields th <C T^p (th is the decoherence time due 
to the magnetic field) in which one can calculate 
the subleading correction to Sa perturbatively in 
the interaction. In this case the dcphasing time 
can be extracted from a perturbative expression 
for 5a{H) only if one makes an additional assump- 
tion about the explicit form of the phase relax- 



ation in time. No such assumption is needed within 

our analysis. In addition to that, the perturbative 
results [11] were presented only in the high tem- 
perature limit Tth '> 1 which, being combined 
with Th '^Tip, is equivalent to a strong inequality 
Tt^ >>> 1. No such inequality was imposed in 
our calculation [5] . 

In this paper we establish a direct relation be- 
tween perturbative and nonperturbative results in 
the problem of quantum dcphasing in disordered 
conductors. We will demonstrate that (i) the per- 
turbation theory in the interaction is insufficient 
for the problem in question and (ii) in order to de- 
termine Tip it is necessary to perform a nonpertur- 
bative analysis with the exponential accuracy. We 
analyze both the dcphasing rate and the system 
magnctoconductance for various temperatures and 
the magnetic fields. 



2. Insufficiency of the perturbation theory 

Originally the dcphasing time of electrons in 
disordered conductors was introduced and eval- 
uated by Altshulcr, Aronov and Khmclnitskii 
within a phenomenological nonperturbative pro- 
cedure [2]. Later Fukuyama and Abrahams [17] 
studied the problem perturbatively in the inter- 
action. Since then the perturbative approach was 
frequently used in the analysis of quantum dcphas- 
ing. We are going to demonstrate that this pertur- 
bative strategy is insufficient and does not allow 
to unambiguously determine the dcphasing time 
even in those limits where the weak localization 
correction to the conductance can be calculated 
perturbatively in the interaction. 

In what follows we will restrict our analysis to 
quasi-ld disordered systems. Extension of our re- 
sults to higher dimensions is straightforward [18]. 
In the Id case the weak localization correction to 
the wire conductance reads 

Sa.iH) = J ^e-*/-F(Vr,). (1) 
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Here the function F{t/T^) accounts for the 
electron-electron interaction {F = 1 without in- 
teraction). Provided in the long time limit the 
function F decays faster than the integral 

(1) converges even for 1/th = and we get 



6ai = -a—jDT^, (2) 

TT 

where a ^ 1 depends on the particular form of the 
function F . This result is obviously nonperturba- 
tive in the interaction. E.g. if one formally defines 
the "interaction strength" A in the effective Hamil- 
tonian, in certain limits one will obtain I/t^ oc A 
and, hence, 8a\ (x I/a/A. For tr > Tip any attempt 
to expand the conductance in powers of A is mean- 
ingless and may only yield divergences in all orders. 

In order to avoid this problem the authors [11] 
suggested to consider the limit of strong magnetic 
fields th <C t^. In this case the integral (1) is cut 
at times t ^ th much shorter than and the 
expansion of dai in powers oi tr/tip oc A can be 
performed. Keeping only the leading terms of this 
expansion, from (1) one readily finds 

5ai -da\' ^ F (0), (3) 

-iTT Ttp 

where 5a f'^ = —{e^/-K)yJDTH- We observe that a 
perturbative expansion in the interaction allows us 
to determine the combination F'{fS)/T^ but not t<^ 
itself. The full function F{t/T^) (and thus F'(0)) 
will remain unknown. Although in the limit th <C 
Tip the value 5a\{H) can be calculated perturba- 
tively in the interaction, this would yield no in- 
formation about the dephasing time t<^. Such in- 
formation can be extracted from the perturbation 
theory only if one assumes some particular form 
of the function F{t/T^). But this form should be 
found rather than assumed. This task can be ac- 
complished only if one goes beyond the perturba- 
tion theory. 

In order to illustrate this conclusion let us con- 
sider several different functions F{t/T^). Perhaps 
the most frequent choice of this function is 

Fi(f/r^) =exp(-VT^). (4) 



This assumption was also adopted in [11]. Perform- 
ing a perturbative calculation in the high temper- 
ature limit Tth ^ 1 the authors [11] arrived at 
the result Sai — Sa'f^ oc Tr^ which, being com- 
bined with (3) and (4), yields the dephasing rate 

[11] 1/Tpl OC T^/t^. 

Another possible choice of the function F can be 

F^{t/T^)=eM-it/Tpf'^). (5) 

For Th > Tp this function - as well as (4) - yields 
eq. (2). However, in the limit th ^ T,p from (1) 
and (5) one finds 

Clearly, this result is incompatible with eq. (3) be- 
cause F'{0) should not depend on th and Tp. At 
the same time from (5) one has F'{0) = 0, and 
again the contradiction between (3) and (6) is ob- 
vious. Also, combining (6) with the perturbative 
results [11] one arrives at 1/Ty2 « T^/s agree- 
ment with [2,3] but in a clear disagreement with 
1/t,^i found in [11]. 

Finally, let us choose the trial function F in the 
following form: 




where 6^1. For th > Tp, the result (2) is recovered 
again, while in the limit th Tp one arrives at 
eq. (3) with F'{Q) — 6/4—1. Combining this result 
with those of [11] one finds Tip3 oc (4 — b)/{TT]j'^). 
For 6 < 4, & = 4 and h > A this result would im- 
ply respectively positive, infinite and even negative 
dephasing rates, which is an obvious nonsense. 

The above examples clearly demonstrate that a 
perturbative procedure is principally insufficient 
for our problem because it leads to completely 
ambiguous results for t^,. At the same time the 
nonperturbativc cq. (1) yields practically indistin- 
guishable magneto conductance curves (see Fig. 1 
of [18]) and the same dephasing time Tp, (up to 
a prefactor a ~ 1 which can be absorbed in Tp, 
anyway) for all the trial functions Fi, F2 and F3. 



3 



In order to obtain correct information about Tip 

one should determine the function F at times t ^ 
Tip. This can be done only by means of a nonpertur- 
bative calculation simply because for t ^ Tip there 
exists no small parameter in the problem. Choos- 
ing the limit tu <C t<^ enables one to find bo\ (H) 
but not Tip. We can also add that the phenomeno- 
logical procedure [2] allows to non-perturbatively 
treat only the contribution of "classical" modes 
cu < T, whereas the effect of "quantum" modes 
a; > T is accounted for by means of the nonpertur- 
bative analysis [5]. This analysis will be extended 
further in the next section. 



3. Nonperturbative results 

Let us express the function F (1) in the form 

Fit)=A{t)eM-m) (8) 

Without interaction one has A{t) = 1 and f{t) = 
0. For a complete description of the effect of in- 
teraction on the weak localization correction (1) it 
is in general necessary to evaluate both functions 
f{t) and A{t). An important observation is, how- 
ever, that information about the effect of interac- 
tion on A{t) is not needed to correctly evaluate the 
dephasing time Tip, it is sufficient to find only the 
function f{t) which describes the decay of correla- 
tions in time and provides an effective cutoff for the 
integral (1) at f ~ Tip. The role of the pre-exponent 
A{t) is merely to establish an exact numerical pref- 
actor. Since Tip is defined up to a numerical pref- 
actor of order one anyway, it is clear that only the 
function f{t) - and not A{t) - is really important. 

3.1. Exponent 

The function f{t) can be straightforwardly eval- 
uated by means of the path integral formalism. 
This procedure amounts to calculating the path in- 
tegral for the kernel of the evolution operator over 
the particle coordinate r and the momentum p [5] 



J ~ j Dr j Dp cxp{iSo - iS'o - iSr - Si) 



(9) 



within the saddle point approximation on pairs of 
time reversed paths and to averaging over diffusive 
trajectories. Here So and S'q represent the electron 
action on the two parts of the Keldysh contour, 
while iSn + Si accounts for the interaction. The 
general expressions for Sr and Si were derived in 
[5] . Within RPA these expressions contain the full 
information about the effect of interaction in all 
orders of the perturbation theory. 

The saddle point approximation procedure was 
also described in details in [5]. One can demon- 
strate that the contribution of the real part Sr of 
the action vanishes on any pair of diffusive time 
reversed paths. Note that such cancellation is a 
generic property of a wide class of influence func- 
tionals describing dissipative environments. E.g. 
the same cancellation is observed in the Caldeira- 
Leggett model [7] and in some other models [9]. 
Hence, the function f(t) in the exponent (8) is de- 
termined only by the imaginary part of the ac- 
tion Si [5] . For the sake of generality here we will 
present the resiilt for all dimensions. Evaluating 
Si on pairs of time reversed saddle point paths we 
find [18] 

4g2£,l-d/2 f ^d^ f doj dio' {1- COS LOt) 



^^^>~ ad{2T,Y J 1 + x^J 



.'\d/2-2 



{OJ-LU') COth ^ 



2T 



+ 



1I/2-2 



+ - 



' COth ■ 



2T 



• (10) 



Note that in Id and 2d cases the integral of the 
first term over u' diverges at uj' 0. However, this 
divergence is exactly canceled by the second term 
and the whole integral is finite in any dimension. 

Evaluating the result (10) for a Id case in the 
quantum regime irTt <C 1 we obtain 



2D 2e^ Dt / 27rt \ , , 
/ \ —t+ \/— In 6 ,11 
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where Tg = Wir/Zis the elastic mean free time. Note, 
that apart from the leading linear in time term 
there exists a smaller term oc \/tln{t/Te), which 
also grows in time. In the opposite thermal limit 
nTt > 1 eq. (10) yields 




V2 



(12) 



where C,{x) is the dzeta-function. We observe that 
in both cases (11) and (12) there exists a linear 
in time temperature independent contribution to 
f{t) which determines the dcphasing time at low 
temperatures [5,6]. Beside that a.t Tt ^ 1 there 
exists another term oc Tt^^"^ which yields dominat- 
ing contribution to at high temperatures T > 
Tq ^ l^Jiyfl, where the result [2] T^p oc T"^/^ is 
recovered. 



3.2. Pre-exponent 



and no compensation of the exponential decay oc 
exp(— /(t)) by the pre-exponential function A{t) 
can occur in the long time limit. Another impor- 
tant property [18] is that in the interesting limit of 
low temperatures the effect of interaction on the 
pre-exponent becomes important at i ~ r<^, i.e. 
on the same time scale as for the function f{t) in 
the exponent. These two properties allow to com- 
pletely ignore the effect of interaction on the pre- 
exponent [5]. 

In the short time limit t <^ the correction to 
the pre-exponent due to interaction is small and 
one can proceed perturbativcly in the interaction. 
At T in the leading approximation one finds 
[18] 



A{t) ~ 1 + 




(13) 



This equation is important for deriving the pertur- 
bative expressions for the magnctoconductance in 
the limit t// <C r,^. The corresponding results are 
presented below. 



A rigorous calculation of the pre-exponential 
function A{t) = A{t)/^ in (1) in the presence of 
interactions is beyond the frames of the present 
paper. Fortunately the precise form of A{t) for all 
times is not interesting for us here. Of importance 
is to qualitatively understand how the function 
A{t) is modified in the presence of the electron- 
electron interaction. 

It is well known [3,19] that without interactions 
the function A{t) is related to the return probabil- 
ity of diffusive trajectories to the same point after 
the time t. In the presence of dissipation (described 
by the term Sr in the effective action) the particle 
energy decreases and its diffusion slows down. This 
implies that at any given time t the function A{t) 
should exceed l/Vi- On the other hand, at least 
if the interaction is sufficiently weak, diffusion will 
take place at all times and, hence, A{t) will always 
decay as l/-\/t or slower. 

The latter - intuitively obvious - property of 
the pre-exponent was confirmed by our analysis 
[18]. The function A{t) indeed decays at all times 



4. Perturbation theory for the conductance 

Let us evaluate the weak localization correction 

(I) in the limit of strong magnetic fields th 
Ttp. Performing the short time expansion of both 
exp(— /(f)) and A{t) to the first order and combin- 
ing (1), (8) with (11), (13), in the limit Tth < 1 
we find 

,,,_,,(0)^f!£!?^LfM_3.74V(14) 

It is easy to observe that this result does not con- 
tain the zero temperature dcphasing time at all. 
The linear in time terms in the expressions for f(t) 

(II) and A{t) (13) cancel each other exactly in the 
first order. The same cancellation occurs in the 
hmit Tth > 1- Again combining (1), (8) with (12), 
(13) and expanding exp(— /(i)) to the first order 
in /, at Tth » 1 we obtain 

^a,-^ar^--^Tri[±+ f^'^ 1 .(15) 
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In order to establish the exact numerical prefactor 
in front of the last term in (15) we also took into 
account the T-dependent linear in time contribu- 
tion to A{t) [18] omitted in eq. (13). We note that 
the result (15) agrees (up to a numerical prefactor 
of order one in front of the term cx Trfj) with the 
perturbative result [11] obtained in the same limit 
Tth ^ 1- However, as it was already discussed 
above, no information about the dephasing time at 
low T can be extracted from (15). 

Finally, let us derive the expression for the weak 
localization correction in the limit of high temper- 
atures and weak magnetic fields th ^ t^. In this 
limit the phase relaxation is determined by the 
term fd oc Tt^/"^ in (12). Keeping this term in the 
exponent, expanding the exponent in f{t) — fci{t) 
to the first order and making use of a short time 
expansion of A{t) we get 



(5cti 



1/3 



r(i/3) 



C(l/2) eVD 

(2^)3/2 



.(16) 



Here r(a;) is the Euler's gamma-function. The first 
term in the right-hand side of this expression cor- 
responds to the classical result [2] , while the second 
term represents the first correction to this result. 
With decreasing temperature this correction grows 
faster than the absolute value of the first term. At 
the same time in the high temperature limit the 
second term in (16) always remains much smaller 
than the first one and the higher order terms should 
also be taken into account already at T > Tq. 

Note that the result (16) - as well as eqs. (14), 
(15) - does not depend on the zero temperature 
value of Tip , this value drops out due to the same 
cancellation of the linear in time T-independent 
terms from the exponent and the pre-exponcnt. 
Thus also the first order high temperature expan- 
sion of the weak localization correction to the con- 
ductance cannot provide correct information about 
the electron dephasing time at low temperatures. 

In conclusion, we established an explicit rela- 
tion between perturbative and nonperturbative re- 
sults in the problem of quantum dephasing in dis- 
ordered conductors. The dephasing time remains 



finite down to T = due to the electron-electron 
interactions. 
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